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Abstract 



We study how well KamLAND, the first terrestrial neutrino experiment 
capable of addressing the solar neutrino problem, will perform in as- 
certaining whether or not the large mixing angle MSW solution (with 
Amfi < 10 4 eV 2 and oscillation amplitude sin 2 2 6>i 2 > 0.3) is cor- 
rect. We find that in a year of operation KamLAND will provide unequivo- 
cal evidence for or against this solution. Furthermore, its sensitivity to the 
three-neutrino oscillation parameters in this region is sufficiently acute as 
to determine Am^ to approximately ±10% (for sin 2 2#i2 > 0.7) and to fix 
sin 2 2 6»i2 to within ±0.1 (at the 2a level) with three years of accumulated 
data, independent of the value of #13. 
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I. INTRODUCTION 



Several neutrino oscillation experiments now indicate that neutrinos are massive and 
that neutrino flavor mixing occurs. Atmospheric neutrino experiments (Kamiokande [|TJ, 
SuperKamiokande @, 1MB ||, Soudan [|] and MACRO 0) report a v^/v e event ratio 
that is about 0.6 times the expected ratio; the flux shows a zenith angle dependence 
that is explained by oscillations. All experiments that measure the solar neutrino flux 
(Homestake §, SAGE 0, GALLEX [§, Kamiokande |§ and SuperKamiokande 0) find a 
deficit of 1/2 to 1/3 of the Standard Solar Model prediction ||11|| . LSND JT2|, an accelerator 
experiment, finds evidence for — > v e and — > v e oscillations, a result that is not excluded 
by KARMEN |13| and awaits confirmation by MiniBooNE [Q. The current generation of 
nuclear reactor experiments (Palo Verde |TJ| and CHOOZ [|T^]) find null oscillation results 
and rule out v e — > v x oscillations for Am^ > 10~ 3 eV 2 at maximal mixing and sin 2 2 6*13 > 0.1 
for larger Am^ (at the 95% confidence level). 

The focus of the present work is the solar neutrino puzzle |I7J and KamLAND's |18| role 
in resolving whether or not the large mixing angle (LMA) solution (see e.g. Ref. |19|) is 
the correct one. The LMA solution has Am^ in the range 10~ 5 to 10~ 4 eV 2 and mixing 
amplitude sin 2 2#i2 > 0.3. Data from SuperKamiokande |20| now favors the LMA solu- 
tion over the small mixing angle, low and vacuum oscillation solutions (see Ref. for a 



discussion of the various oscillation solutions). Solar neutrino measurements at SNO f22 



should discriminate between different oscillation solutions, but perhaps not in its first year 



of operation p3| . If the LMA solution is the correct one, KamLAND will provide a precise 
determination of the oscillation parameters. 

KamLAND is unique in its potential as the first terrestrial experiment to explore the 
solar neutrino anomaly. It will provide a definitive test of the LMA solution by either ruling 
it out or by pinning down the values of Am^ and sin 2 2# 12 . KamLAND should be able to 
provide this answer in a year from the start of running in spring 2001 and then provide an 
accurate determination of the solar neutrino oscillation parameters by the end of the three 
years over which it is expected to take data for. Particularly significant will be its ability 
to precisely determine Am^ in the LMA region because unlike the case of solar neutrinos, 
the Am^-dependent contribution to oscillations of the reactor neutrinos will not suffer from 
averaging of the oscillation L/E dependence. 

If the LMA solution is correct, Am^-dependent CP-violating effects [23] can be large 



enough to be tested at very long baseline accelerator-based experiments | 25|| . With Kam- 
LAND promising a precisely known value of Am^, it will provide essential information for 
future experiments studying CP-violation in the neutrino sector. 

In Section II we provide a brief overview of the KamLAND experiment and the oscillation 
hypothesis in which we will work. Section III will be devoted to details of our simulation of 
the experiment and the subsequent data analysis. We conclude in Section IV. 



II. THE KAMLAND EXPERIMENT 

KamLAND is a reactor neutrino experiment with its detector located at the Kamiokande 
site. With low energy neutrinos ((E v ) ~ 3 MeV), it can only measure v e disappearance and 
therefore will be unable to access small mixing angles sin 2 2# 12 < 0.1. About 95% of the v e 
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flux incident at KamLAND will be from reactors situated between 80 — 350 km from the 
detector, making the baseline long enough to provide a sensitive probe of the LMA solution 
of the solar neutrino problem. Specifically, the sensitivity to the mass-squared differences 
will lie between Am 2 2l ~ (L (m)/E (MeV)) -1 > 10~ 5 eV 2 and < 10~ 4 eV 2 . Despite the 
absence of a single baseline, the measured positron energy spectrum allows a sensitive probe 
of oscillation effects. 

We consider a framework of three-flavor oscillations among the massive neutrinos 
(yi, V-1-, Pz). In the phenomenologically interesting limit, Am 2 = m\ — m\ <C m 2 — m 2 ,, 
the survival probability at the detector is given by 

99 994 9 /1.27Am 2 (eV 2 )L(m)\ . , 

P(u e - i/ e ) = 1 - 2 4 c 2 3 - 4 s\ 2 c\ 2 4 sin 2 { e^AbV) ) > W 

where s?- = 1 — c?- = sin 2 ^ are the neutrino mixing parameters, L is the distance of the 
detector from the source and E v is the energy of the anti-neutrino. Here we have averaged 
over the leading oscillations (sin 2 (1.27 \m\ — m\\ L/E u )) = 1/2. Matter effects on oscillations 
are negligible at the KamLAND baseline. 

Although we restrict our analysis to three neutrino oscillations, it is possible that the 
transition v e — > u s , where v s is a sterile neutrino, is responsible for the depletion of the v e 
flux. (See Refs. f26|,|27j for models and analyses of four-neutrino oscillations for the solar and 
atmospheric anomalies). For example, in a 2 + 2 neutrino mixing scheme, in which one pair 
of nearly degenerate mass eigenstates has maximal v e — > v s mixing for solar neutrinos and 
the other pair has nearly maximal mixing — > v T for atmospheric neutrinos, the survival 
probability is p6| 

v . o 2 /1.27Am 2 (eV 2 )L(m)\ 

P(v e -> u e ) = 1 - 4 e 2 - sin 2 ^ { K ' . (2) 

V ; V ^(MeV) / ' y ! 

where e ~ (0.016 eV 2 /Am 2 5JVD ) - 91 is restricted by KARMEN p| and BUGEY @ to 



be in the interval (0.01,0.1) corresponding to a Am 2 LSND range 0.2 to 1.7 eV 2 H). In 
Eq. (Q) we have averaged over the leading oscillations associated with the LSND mass scale. 
On comparing Eq. ([I]) (fixing 6> 12 = vr/4 and sin 2 2 6 , 13 = 0.1) with Eq. (g), one sees that 
KamLAND will be unable to distinguish whether the oscillations are to flavor or to sterile 
neutrinos, but neutral current measurements at SNO will accomplish this. KamLAND is 
more sensitive to #13 than to e, and we find that it will not provide useful information about 
9i 3 and consequently even less about e. In the following, we base our analysis on Eq. ([!]). 

The target for the v e flux consists of a spherical transparent balloon filled with 1000 tons 
of non-doped liquid scintillator. The anti-neutrinos are detected via the inverse neutron 
/3-decay 

v e + p — > e + + n . (3) 
The cross section for this process is |29[ 



a{E v ) = -^—{E v - AM) \{E V - AM) 2 - m 2 



eJ 



0.952 (^-AM)[(^-AM)'-^] x ^ ^ (4) 
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Reactor 


Distance (km) 


Percent of total flux 


Kashiwazaki 


160.0 


33.36 


Ohi 


179.5 


14.76 


Takahama 


190.6 


9.76 


Shiga 


80.6 


8.51 


Tsuruga 


138.6 


8.12 


Mihama 


145.4 


8.10 


Hamaoka 


214.0 


8.06 


Fukushima — 1 


344.0 


4.17 


Fukushima — 2 


344.0 


3.86 


Tokai — II 


294.6 


1.30 



TABLE I. The expected relative contribution of each reactor to the v e flux detected at Kam- 
LAND in the case of no oscillations. Only reactors in a 350 km radius of the detector are considered. 



where m e is the mass of the electron, AM = m n — m p is the neutron-proton mass difference, 
T n = (886.7 ± 1.9) s is the neutron lifetime and / = 1.7152 is the phase space factor which 
includes Coulomb, weak magnetism, recoil and outer radiative corrections. The theoretical 
error in the above cross section is less than a percent. The neutrino capture process has a 
threshold of AM + m e = 1.804 MeV and the v e flux above this threshold is 1.3 x 10 6 cm" 2 s _1 



which is known to a precision of about 1.4% [30|. In the case of no oscillations, KamLAND 



expects to see ~ 800 events per year with a background ~ 40 events per year. The dis- 
tribution of the background events versus positron energy is expected to be known, thus 
facilitating a clean extraction of the signal. We will assume that a background subtraction 
can be made for our data simulation and analysis. 



III. DATA SIMULATION AND ANALYSIS 

To calculate the v e flux we include contributions from all nuclear reactors within a radius 
of 350 km from the detector. Table | gives the relative fluxes (without oscillations) and thus 
shows the relative importance of each reactor to the experiment (total fluxes from each 
reactor are tabulated in Ref. |18|). The v e spectrum above the 1.8 MeV threshold for the 
process ([|) is the result of the decay of fission fragments of the isotopes 235 U, 239 Pu, 238 U 
and 241 Pu. The spectrum from the fission products of each of these isotopes can be found 
3~T|. As the reactor operates, the concentration of 235 U decreases and that of 239 Pu and 



m 

241 



Pu increases. We do not account for the fissile isotope evolution and instead assume a 
typical fraction of fissions (as in Ref. |ftj) from the four fissile materials (see Table |TJ). We 
have confirmed that the effect of the evolution is small by reproducing the v e spectrum from 



the Palo Verde experiment [JT5 



With the relative flux from each reactor known, and the knowledge of the v e spectrum 
(assumed to be the same for all reactors), we calculated the e + spectrum at KamLAND 
resulting from inverse /3-decay and normalized it to yield a total of 800 events per year in 
the absence of oscillations. This procedure effectively accounts for the effective number of 
free protons in the target, the e + and n efficiencies, the efficiency of the e + — n distance cut 
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Start of Cycle (%) End of Cycle (%) Typical value (%) 



235 jj 


60.5 


45.0 


53.8 


239p u 


27.2 


38.8 


32.8 


238 jj 


7.7 


8.3 


7.8 


241 Pu 


4.6 


7.9 


5.6 



TABLE II. The fraction of fissions from the four fissile elements in a nuclear reactor at the 



beginning of the cycle, the end of the cycle and a typical value during a cycle taken from Ref. [ 30 ] . 



and fluctuations in the power output of the reactors arising from dead time for maintenance 
and seasonal variations of power requirements. 

Figure [I] shows the e + energy spectrum expected at KamLAND with three 
years of data, illustrated for the cases of no oscillations (dotted histogram) and 
(Am 2 , sin 2 2 #i2) = (7 x 10 -5 eV 2 , 0.75) (solid histogram). Each simulated data point is gener- 
ated by randomly choosing a point from a Gaussian distribution centered at the theoretical 
value and of width equal to the square root of the theoretical value. The lower plot shows 
the ratio of the simulated data to the expectation for no oscillations. The errors shown 
are statistical. The plots are versus the total e + energy. A plot of the total measurable 
energy would be shifted to the right by m e = 0.511 MeV because the v e signature involves a 
measurement of the kinetic energy of the e + and the annihilation energy in the form of two 
0.511 MeV gamma rays. Figure |2] illustrates the significant changes in the e + spectrum for 
different values of Am 2 at a given value of sin 2 2 6 l i2. The spectra overlap to a large extent 
for values of Am 2 > 2 x 10" 4 eV 2 . 

For the statistical analysis we define x 2 as 

2/ A 2 • 2r> n \ V» (Simulated ~ A^, Am 2 , shl 2 2 12 )) 2 

X (Am , sin 2 012) = 2^ ttt, , (5) 

i = l ! 'simulated y^i ) 

where % labels the 17 e + energy bins each of width 0.4 MeV and midpoint A r simmat ed(-E'0 is 
the number of simulated events in the i th bin and N(Ei, Am 2 , sin 2 2 Oyz) is the corresponding 
theoretical value for oscillation parameters (Am 2 , sin 2 2 612). The only fitted parameters are 
Am 2 and sin 2 2#i2. We keep sin 2 2#i3 fixed when performing the \ 2 analysis (using the 
Levenberg-Marquadt method). From CHOOZ we know that sin 2 2# 13 < 0.1 and we consider 
the two extreme cases sin 2 6 13 = and sin 2 2# 13 = 0.1. 

It is important to note that the main source of uncertainty in the experiment comes from 
conversion of the fission rates in the reactors to v e fluxes and the normalization uncertainty 
is expected to be less than 3% |32|. (The CHOOZ experiment [16| had a normalization 



uncertainty of 2.7%). Even if the combined systematic uncertainty could be as large as 5%, 
the shapes and sizes of the confidence contours in our analysis would be unaffected by its 
inclusion. We can thus safely ignore systematic uncertainties in what follows. 

We concentrate on the region defined by sin 2 2 12 > 0.2 and 10" 5 < Am 2 < 2 x 10" 4 eV 2 . 
For values of Am 2 outside this region, sensitivity to the neutrino energy-dependence is lost. 
For Am 2 < 10~ 5 eV 2 , the value of the oscillation-dependent sinusoidal factor gets small. 
For Am 2 > 2 x 10~ 4 eV 2 , the argument of the Am 2 -dependent sine function in Eq. ([]]) 
becomes large because of the long baseline and the oscillations get averaged. By studying 
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the behavior of y 2 — Xmin as a function Am 2 for several simulated datasets, we find that for 
theoretical inputs with Amf heor < 2 x 10~ 4 eV 2 , x 2 ~ Xmin h &s unique, well-defined minima 
at the theoretical values. For inputs with Amf heor > 2 x 10 _4 eV 2 we find the presence of 
more than one acceptable minimum in the region 2 x 10~ 4 < Am 2 < 10~ 3 eV 2 and a 
continuum of solutions Am 2 > 1CT 3 eV 2 that are acceptable at the 2a level. Thus, for values 
larger than 2 x 10 _4 eV 2 , KamLAND can place a lower bound on Am 2 , but cannot uniquely 
determine the value. However, the region of good sensitivity covers the entire LMA solution. 

Figure |3| shows fits to e + spectra (expected after three years of running) for values of 
Am 2 and sin 2 2 6 l 12 covering the entire region of the LMA solution with 9 13 = 0. The plot 
shows lcr (68.3%) and 2a (95.4%) confidence contours. The diamond is the theoretical value 
for which data was simulated and the cross marks the best fit point. Each point is labelled 
by the expected number of signal events. The x 2 value for the best fit corresponds to a x 2 
probability of at least 80%. For a given value of Am 2 , the confidence regions get flatter along 
the Am 2 -direction as sin 2 2# 12 is increased from 0.2 to 1, resulting in a finer determination 
of Am 2 for larger sin 2 2 # 12 . The value of sin 2 2 # 12 will be determined to an accuracy of ±0.1 
(2a) throughout almost all of the parameter space under consideration. The Am 2 value will 
be determined within ±10% (2a) of its actual value for 10 -5 < Am 2 < 2 x 10~ 4 eV 2 and 
sin 2 2fl 12 > 0.7. 

We found that for sm 2 2 9 13 = 0.1 the contours are almost identical to that of Fig. |3|; 
thus, KamLAND will not be sensitive to #13 in any region of the parameter space. In Fig. |] 
we show the expectation with one year of accumulated data (taking #13 = 0) and find that 
in this time period KamLAND should already give us the LMA parameters to reasonably 
good accuracy. 



IV. CONCLUSION 

The results of our work are succinctly summarized in Figs. Within one year of 

operation, KamLAND measurements of the reactor v e flux will establish whether or not the 
LMA solution (10 -5 < Am 21 < 10~ 4 eV 2 and sin 2 2 6> 12 > 0.3) is correct. If it is, in a span 
of three years the accuracy with which sin 2 2# 12 , and especially Am 2 1; will be determined 
is striking. KamLAND will give us sin 2 2 6 l i 2 to an accuracy of ±0.1 and Am 21 to within a 
factor of 2 for sin 2 2 9u > 0.2; Am 21 will be known to an accuracy of ±10% for sin 2 2 Qyi > 0.7 
(at the 2a level). However, no knowledge will be gleaned about #13. Our conclusions assume 
that the experiment employs a 1000 ton fiducial volume. If the fiducial volume is instead 
600 tons, Fig. f| applies for data collected over 20 months instead of 12 and the contours 
of Fig. [3] will be roughly twice as large at the end of three years. Details of systematic 
uncertainties and the accuracy with which the background shape is known will cause minor 
changes in our results, but the overall conclusions will remain intact. 

If the LMA solution is found to be correct, the precision with which Am 21 will be pinned 
down at KamLAND will prove to be very beneficial to studies of CP- violation in the lepton 



sector in long baseline experiments [pi, 35 
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sirL 2<9 12 = 0.75 



LhUT 
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e + Spectrum 



Simulated Data With. Oscillations 
Expectation Without Oscillations 



4 



Total e + Energy (MeV) 
FIG. 1. The expected e + energy spectrum (versus total e + energy) with three 
years of accumulated data for the case of no oscillations (dotted histogram) and for 
(Am 2 , sin 2 2 6u) = {7 x 10 -5 eV 2 , 0.75) (solid histogram). The data points represent the simulated 
spectrum. The lower plot shows the ratio of the simulated data to the expectation for no oscilla- 
tions. The errors shown are statistical. 
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Total e + Energy (MeV) 
FIG. 2. KamLAND's sensitivity to Am 2 is unprecedented. In three years it will easily be 
able to discriminate between only slightly different values of Am 2 in the LMA region. For 
Am 2 ^ 2 x 10~ 4 eV 2 the spectra overlap significantly. 
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FIG. 3. Fits to e + spectra for values of Am 2 and sin 2 2#i2 covering the entire region of the 
LMA solution with #13 = 0. The la (68.3%) and 2a (95.4%) confidence contours are shown. 
The diamond is the theoretical value for which data was simulated and the cross marks the best 
fit point. Each point is labelled by the expected number of signal events in three years. If no 
oscillations occur, the expectation is 2400 events. 
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FIG. 4. The same as Fig. [3| except that only one year of accumulated data is assumed. Cor- 
respondingly, each point is labelled by the expected number of signal events in one year. If no 
oscillations occur, the expectation is 800 events. 
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